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Abstract 

We propose a scheme for symmetrization verification in two-particle 
systems, based on one-particle detection and state determination. In con- 
trast to previous proposals, it does not follow a Hong-Ou-Mandel-type 
approach. Moreover, the technique can be used to generate superposition 
states of single particles. 

PACS: 03.65.-w; 03.67.Mn 

1 Introduction 

Non-factorizable states play a fundamental role in quantum theory. There are 
two classes of non-separable states, entangled and (anti)symmetrized ones. Re- 
cently, there has been a boom in the theoretical and experimental study of the 
first ones, specially concerning their generation, verification and potential ap- 
plications in the fields of quantum information and conceptual foundations of 
quantum theory. Here, we are mainly interested into the question of verification. 
References [1, 2] are two good introductions to the problem. 

In contrast, there are not so many studies concerning the same issue for 
(anti)symmetrized states of massive particles. In quantum theory, two identical 
particles whose wave functions overlap at a given instant are indistinguishable at 
any subsequent time. If one localizes one of the particles at any instant after the 
overlapping, it is impossible to identify it with any of the initial ones. In order 
to correctly describe this impossibility the wave function of the two-particle 
system cannot be a product state (even for non-interacting particles) but must 
be symmetrized in the case of bosons and antisymmetrized for fermions. 

In quantum optics, the standard test for determining identicalness, i. e., for 
testing if two photons are really in indistinguishable states (and consequently it 
is necessary to symmetrize the state) is based on the Hong-Ou-Mandel (HOM) 
arrangement [3] . This technique is based on the interference of the photons at 
a beam splitter [4]. 

Several authors have proposed to use the same approach for massive particles 
[5, 6, 7, 8]. In [5, 6], mesoscopic electron beam splitters based on electron beam 
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lithography on a GaAs high-mobility two-dimensional electron gas system have 
been experimentally used to test the antibunching of fermions, an unequivocal 
demonstration of their antisymmctrization. Several other arrangements can be 
used as beam splitters for beams of massive particles [9, 10], providing also 
an alternative basis for feasible massive HOM-type tests. Unfortunately, their 
manipulation is much more delicate than their counterparts for light beams, 
making more difficult the application of the technique to massive particles. 

The importance of having methods able to verify (anti)symmctrization can 
be easily understood. For instance, in [11] it has been suggested to study mas- 
sive two-particle intcrfcrometry by diffraction gratings. One expects novel ef- 
fects when identical particles arrive in (anti)symmetrized states. In order to 
experimentally test such effects we must know in advance if the particles in- 
cident on the diffraction grating are in (anti)symmetrized states. If not, a 
negative result in the experiment could be imputed either to the absence of 
the effect or to the state of the particles (being in a product state instead of a 
(anti)symmetrized one). In contrast, if we know by sure that the particles are 
in a (anti)symmetrized state a negative result could only be interpreted as a 
demonstration of the absence of the effect. 

We propose here a scheme able to verify (anti)symmetrization in two-particle 
systems, only based on the detection of one of the two particles and the deter- 
mination of the state of the surviving one. To be concrete, if the state of the 
two-particle system is (anti)symmctrized the mode distribution of the surviving 
particle will be a combination of the initial distributions of the two particles. 
The scheme is, in contrast to previous proposals, a test of identicalness not 
based on a HOM-type approach. 

As an interesting byproduct, the technique also provides an efficient method 
to generate one-particle superposition states. The surviving particle is in a 
superposition of the initial states of the two particles. The technique is specially 
well-suited to generate superpositions of multi-mode states peaked around very 
different momentum values. 



2 The arrangement 

We consider the same arrangement as that in [12], where the process of detec- 
tion of one of the particles in a system of two particles in multimode states is 
analyzed, with special emphasis on the one-particle interferences present in the 
detection pattern. 

Two identical particles are prepared in the state 

\2 fg >= J d 3 p J d 3 q/(p) 3 (q)a+a+|0 > (1) 

where / and g are the mode distributions, which we assume to be normalized, 
/ d 3 p\f(p)\ 2 = 1 = / d 3 p\g(p)\ 2 . In Eq. (1) a+ is the creation operator of the 
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mode p and |0 > is the vacuum state. The creation and annihilation operators 
obey the usual relations [d p , a+] T = <5 3 (p — q). In the double sign expressions 
the upper one refers to bosons and the lower one to fermions. Note that the 
above state is not normalized. We shall work with a non-normalized state (see 
in [12] how this state is normalized), postponing the normalization to the final 
stage of the calculation. 

In all the above expressions we have not included the spin indexes. In order 
to simplify the presentation we assume the two particles to be in the same spin 
state making superfluous the presence of the indexes. In terms of the quantum 
description this assumption is equivalent to have for the spin part of the state the 
form \s >i \s >2, that is, a product state with 1 and 2 labeling the two particles. 
We could also have (for both bosons and fermions) a symmetric state for the 
spin variables, (|s >i \s' >2 +\s' >i \s >i)l\[2. With any of the two above 
choices the part of the state related to the spatial variables must be symmetric 
for bosons and antisymmetric for fermions. We do our presentation in this 
framework, symmetrized states for the spatial or momentum wave functions of 
bosons, and antisymmetrized for fermions. The extension to the case of particles 
in antisymmetrized spin states is straightforward, as it is only necessary to 
consider antisymmetric and symmetric spatial states for, respectively, bosons 
and fermions. 

The two-particle system now interacts with a detector. We only consider 
the cases in which only one detection occurs. Moreover, we assume that in 
the detection process the particle is removed, for instance, by absorption. In 
order to postselect these cases one must use a detector able to distinguish events 
with one- or two-detection events (see [13, 14] for an optical detector with this 
property). 

Next, we show that the detection process is described by the action of the 
Schrodinger field operator. From a mathematical point of view this operator 
t/j(R,t) represents the annihilation of a particle at point R and time t [15, 
16]. If the initial state of the system is a n-particle state \n >, the action 
of the operator leaves a n — 1 particle one, i/>|n >^ \n — 1 >. Thus, this 
operator represents in a natural way processes in which particles are removed. 
Intuitively, one can identify this removal process with the destructive detection 
assumed above. With this identification V>(R, t)\2f g > describes the state of 
the resultant one-particle system. At this point a criticism to this identification 
could be made. Being the detection process a non-unitary one, perhaps we 
should use a mixture instead of a pure state to describe the final one-particle 
system. We argue here favoring the pure state choice. Our argument is based 
on the analogy to detection through photomultipliers in quantum optics [4] (see 
also [17] for a detailed treatment of absorption in multiparticle massive beams 
and its description in the framework of the Glauber scheme). In this type 
of destructive detection the photons are absorbed by the atoms composing the 
photomultiplicr. The description of the interaction between detection device and 
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radiation is done in terms of the positive frequency part of the electromagnetic 
field operator E+. If we denote by \R > the initial state of the radiation field, 
its final state is assumed to be E + |i? >, a pure state not a mixture. From 
that pure state one can evaluate the probability of transition to a state \Rp >, 
as | < Rp\E + \R > | 2 . By analogy with the extensively verified scheme of 
detection by photomultipliers we assume that after the destructive detection we 
have the pure state ^>(R, t)\2f g >. As we shall see in the following sections, this 
assumption gives rise to some conclusions that can be experimentally tested, 
making the assumption verifiable. 

In the expression of the field operator R refers to the point of detection. 
We use this notation to emphasize that we are considering the detector to be 
at a fixed position in each repetition of the experiment. We want also remark 
that it does not act as a physical variable (as it is the case of r) but as a 
parameter. This is a necessary consistency condition for the calculation, because 
in the second quantization formalism there cannot be explicit dependence on 
the spatial variables, unless it is in a parametric form [15]. 

The field operator is given by tp(R,t) = J d 3 ptp p (IL)a p (t) 7 with -0 P a com- 
plete set of orthonormal stationary wavefunctions. The most common choice 
for this set is that of plane waves Vp(R) = (27r?i)~ 3 / 2 exp(ip • R/h). The time 
dependence is carried by the annihilation operator. For planes waves it is given 
by o, p (t) = a p exp(—iE p t/h) with a p the annihilation operator at time t = 
and E p = p 2 /2m the energy. 

From now on, by the matter of simplicity, we shall restrict our considerations 
to stationary situations. Consequently, we can drop the time variables from all 
the expressions. 

3 Mode distribution of the surviving particle 

In our case, the final state after one detection is 

Vi(R)|2 /s > (2) 
The evaluation of the state of the surviving particle is simple 

^(R)|2 /9 >= J d 3 p(^(R) 3 (p)±^(R)/(p))a+|0> 

= V/(R)|l s >±V 9 (R)|l/> (3) 

with V/(R) = Jd 3 p/(p)^ p (R) and \l f >= J d 3 pf(p)a+\0 >. 

The above state has not yet been normalized. We denote by \lh > the 
normalized state of the surviving particle, which is given by \lh >— V>(R)|2/ S > 
/N, with 

7V 2 (R) = |^(R)| 2 + |^(R)| 2 ± 2fie(^(R)V s (R) < l ff | 1/ >) (4) 
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Finally, the mode distribution of the surviving particle is 



The mode distribution is neither / nor g, but a combination of them. The 
coefficients in the combination are given by some functions depending on the 
values of the initial wavefunctions at the point of detection. 

Note that only in the cases where ipf(R) — or ^ g (R) = 0, the surviving 
particle has a mode distribution similar to one of the initial ones \l g > or 
|1/ >. This behaviour can be easily understood. Let us consider one of these 
nodal points, for instance, tp g (R) = 0. Then we know by sure that the particle 
detected is of the type / and, consequently, that the surviving one is of type g 
(h(p)=g(p)). 

In the points where simultaneously we have ipf(R) = and ip g (R) = the 
expression for > is undefined, reflecting the property that at these points 
the probability of detection of the particles is null and it does not make sense 
to speak about the surviving one. 

The fact that the mode distribution of the surviving particle becomes a com- 
bination of the initial ones is a clear manifestation of a superposition process. 
This follows directly from the last term in Eq. (3), but can easily be visu- 
alized if we consider the first-quantization version of the problem, where the 
wavefunction is given by 

^(r) = J d*ph(p) cxp(*p • r/h) = ^y^(r) ± Jf^M*) (6) 

showing the superposition of ?pf and tp g . The wavefunction of the surviving par- 
ticle is a superposition of the two initial wavefunctions. This superposition can 
be physically understood in terms of the indistinguishability of the alternatives 
available to the system, such as stated by the general criterion for quantum su- 
perpositions. This criterion states that when the evolution of a quantum system 
is compatible with several indistinguishable alternatives the state representing 
the system is a superposition of them. In our case, outside the nodal points, 
the detector is unable to provide information about the mode distribution of 
the particle measured, making the alternatives detected particle with mode dis- 
tribution f and detected particle with mode distribution g indistinguishable. 

When the particles are not in (anti)symmetrized states the above results do 
not hold. In the second-quantization formalism used in this paper the absence 
of (anti)symmetrization translates into the use of (anti) commuting creation and 
annihilation operators a and b. The initial state (now normalized) is 

|2/ 9 >ns= J d 3 p j d 3 q/(p) 3 (q)a+S+|0 > (7) 
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If we place a detector of particles of type a (the reasoning for type b is similar) 
the state after one-detection event is 

V>( a )(R)|2 /s>Ars =^ a) (R)|l s >W (8) 

The above expression shows that now the mode distribution of the surviving 
particle is identical to the initial one. 

4 The scheme 

The above results provide the basis for a scheme of (anti)symmetrization verifi- 
cation between the two initial particles. The scheme has two steps. In the first 
stage, the two particles interact with the detector. We postselect all the cases 
of one-detection events. Detectors able to distinguish one- and two-detection 
events would greatly simplify the task. As signaled before, in the optical case 
such types of detectors are available [13, 14]. In the second step, we deal with 
the postselected surviving particles. We must test if the particles are in one 
of the two initial states or in a superposition of them. To this end we can use 
standard tomography techniques or study the spatial detection patterns of the 
surviving particles. If the detectors used in the first step arc unable to dis- 
criminate between one- and two-detection events the postselection must also be 
carried out in the second stage. Then we choose the cases in which, in each 
repetition of the process, there is an event in the detector and another in the 
device testing the superposition. As the efficiency of the detectors of massive 
particles is in general high the scheme seems to be, in principle, an effective one. 
These are the general lines of the scheme. Now we discuss in detail a particular 
realization that is, in principle, viable with present technology. We present it in 
Fig. 1. 

We first consider the stage of preparation of the state. We consider an 
arrangement where each particle is emitted by a different source. Two main 
types of sources have been used in the study of (anti)bunching effects in massive 
particles, magnetic traps [18, 19] and optical lattices [20, 21]. Sometimes the 
traps and lattices use previously generated Bose-Einstein condensates (BEC), 
and in other occasions gases produced by evaporative and sympathetic cooling. 
In [19, 20, 21] the atoms are released at the switch-off of the potentials retaining 
the atoms in the trap or lattice. A number of atoms of the order of 10 4 to 10 5 
are released, making the technique useless for our purpose of preparing two- 
particle systems. From our point of view, more interesting is the method of 
[18], where a weak continuous microwave field is applied to the BEC in the 
magnetic trap for output coupling an atom laser. Some atoms inside the trap 
are spin-flipped and do not longer experience the magnetic potential, leaving 
the trap. The number of released atoms is much smaller than in the previous 
techniques. As a matter of fact, in [18] a mean atom number very close to 2 
was reported, making possible the individual detection of the atoms. With this 
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method, adequate values of the microwave fields, and synchronization of the two 
fields acting on different traps it seems possible to create two- particle states. 

For different parameters of the two traps (different cooling,...) we can have 
particles with different momentum distributions. Another way to generate dif- 
ferent momentum distributions is to collimate the released beam via delimiters 
of different apertures. The uncertainty relation tell us that if the width of the 
aperture is I the momentum spread will be Sp ~ Ti/l. 

When the mean separation between the emission of successive pairs of parti- 
cles is large enough the postselection process of pairs can be done in a simple way 
by comparing the detection times in the detector and the tomographic device. 



• 


• 


V J 






Figure 1: Schematic representation of the proposed arrangement. The boxes 
containing the red and blue circles represent the two sources of particles, with 
typical mode distributions / and g. The black half-circle corresponds to the 
main detector. Finally, the green boxes in the bottom of the figure represent 
the detectors array. 

Next, we must perform tomography of the surviving particle in order to 
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determine its state. As it is well-known [22], the first question in quantum to- 
mography is to determine if the system is in a pure state or it must be described 
by a mixture. In the first case, measuring the detection probability in different 
planes along the beam propagation direction we can determine the complete 
(modulus and phase) wave function [22]. In the second case, only the density 
matrix of the system can be obtained from the measurements. Then our first 
task is to determine if we are dealing with a pure state or a mixture. In Sect. 2 
we have assumed that after the detection the two-particle system collapses into 
a pure one-particle state. A possible way to assess if the surviving particle is 
in a pure superposition state consists in the comparison of the spatial detection 
patterns of the surviving and separately emitted particles. To this end we place 
a detectors array downstream of the main detector (that absorbing one of the 
particles). With this array we can determine the spatial detection pattern. For 
the surviving particle in a pure state we have the distribution 

|^(r det )| 2 - MR)| 2 |^(r det )| 2 + K(R)|V/(r det )| 2 

±2Re(a* f (R)^ (r det )a g (R) V / (r de t)) (9) 

with an obvious notation for the a coefficients, and rdet denoting the location 
of the detectors composing the array. 

In a different run of experiments we obtain \tpf(rd et )\ 2 , \ip g (rd e t)\ 2 , |V>/(R)| 2 
and |^ g (R)| 2 by considering situations where only one of the sources emits parti- 
cles. Then we compare \ip h (r de t)\ 2 with \ipf(R)\ 2 \ip g (r de t)\ 2 + \ipg(R)\ 2 \ipf(r det )\ 2 
in order to see if both distributions coincide (there is no superposition of the 
initial states in the surviving particle, corresponding the situation to a mixture) 
or not (there is a superposition associated with a pure state). Note that we do 
not need to determine the normalization factor N(R) because it is not necessary 
to see if both distributions have the same analytical form. 

With the above procedure we can determine the presence of the superposition 
in the surviving particle and, consequently, if it is in a pure state. If it is so, we 
can complete the tomographic process by position measurements of the particle 
in different planes after the main detector (a detectors array at each one of these 
planes). 

It must be noted that if we only want to know if the two-particle state 
was (anti)symmetrized previous to the interaction with the main detector, it is 
enough to carry out the above procedure of comparing the \tph\ 2 and \ipf\ 2 and 
IV'gl 2 distributions, without need of a much more demanding complete process 
of tomography. In effect, when the surviving particle is in a superposition state 
the initial two-particle wave function necessarily had to be in a (anti)symmetric 
state. The situation is different if, in addition to determine the presence of 
exchange processes in the two-particle system previous to the interaction with 
the main detector, we want to know if the spatial part of the state of the system 
was in a symmetric or antisymmetric state. As discussed before, the bosonic 
and fermionic particles can be in a symmetric or antisymmetric spatial state 
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depending on the antisymmetric or symmetric character of the spin part of 
the state. The symmetric or antisymmetric character of the spatial part of 
the two-particle wave function can be determined from the positive or negative 
sign in the superposition of the surviving particle. However, the value of this 
sign cannot be determined with the method of comparison of \tph\ 2 and |^/| 2 
and IV'gl 2 presented in the previous paragraphs. We would need to know the 
complete expression of iph (or equivalently, the knowledge of -0/ and ip g and the 
use of Eq. (9)). We would need a complete tomography process. 

5 Conclusions 

We have presented in this paper a scheme for (anti)symmetrization verification. 
Our proposal is not based on HOM-type arrangements. The ideas presented 
here seem, in principle, accessible to experimental scrutiny. 

The scheme could also be used to generate superposition states in one- 
particle systems, starting from two-particle ones. If the initial two-particle state 
is (anti)symmetrized, the detection of one of the members of the pair outside 
the nodal points leaves the other in a superposition of the two initial states. 
We have generated superposition via (anti)symmetrization and detection. An 
interesting application of the method would be to generate a multi-mode su- 
perposition state of two very different central momenta. We only need (i) to 
prepare the initial particles with momentum distributions f(p) and g(p) peaked 
around two values pf and p g with \pf — p g \ 3> cr, where a is the typical width of 
the distributions, (ii) to get a non-negligible overlapping of the two particles in 
order to have a (anti)symmetrized state, and (iii) to postselect the cases with 
one-detection events. 

Finally, we must consider the possibility of applying the scheme in quantum 
optics. The HOM technique with beam splitters is very efficient for optical sys- 
tems, making superfluous the existence of other verification schemes. Neverthe- 
less, our proposal could be interesting to generate superpositions of one-photon 
states. A simple calculation shows that a result similar to that obtained for the 
distribution of the surviving massive particle holds in quantum optics. More- 
over, in this field there is much experience in dealing with two-photon states, 
even in multimode states [4]. 
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